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1 Introduction 

A necessary part of any one-loop calculation is the loop integration over the undetermined 
four-momentum k in each diagram contributing to the process considered. A loop integra- 
tion will become singular if the 4-dimensional real hypercontour over which the integration 
is performed becomes pinched by two (or more) poles associated with the denominator 
factors in the integrand. Such singularities are known as Landau singularities, and they 
have been studied for some time [1]. 

The denominator of a one-loop integral is equal to the product of propagator denom- 
inators in the associated Feynman diagram, which is independent of the nature of the 
particles in the diagram. Thus, the Landau singularities in a particular Feynman diagram 
are independent of the nature of the particles in it. The behaviour of the integral at a 
singular point can however be affected by the nature of the particles in the diagram, which 
determines the numerator of the loop integral. If the numerator vanishes at the singular 
point, then the integral could be less singular than expected there, or even finite. 

A relevant example of a one-loop calculation in which Landau singularities are encoun- 
tered is gg — > ZZ via massless quark boxes. Three of the six box topologies contributing 
to this process are sketched in figure 1 - the other three only differ by the direction of the 
arrow in the closed quark loop, and give the same contributions as the boxes drawn. 

Apart from mundane threshold singularities, the loops in figure 1 contain Landau 
singularities that are associated with the initial state and loop particles being massless. In 
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Figure 1. Box topologies contributing to gg —> ZZ. 
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fact, all of the diagrams in figure 1 contain at least one of these singularities for arbitrary 
values of the external invariants. Every loop integral contains a collinear singularity (to be 
more precise, two collinear singularities), which is so named because it is associated with a 
quark-antiquark pair attached to one of the external gluons becoming on-shell and collinear 
to that gluon. The first and second loops in figure 1 also contain a soft singularity, which 
is associated with the four-momentum k<i or I2 shrinking to zero. 

We shall refer to the final box topology in figure 1 as the crossed box due to its 
appearance when drawn with initial states on the left and final states on the right (although, 
for reasons of clarity, we will not draw it in this way elsewhere in this paper - see figure 
3 for example). This contains a singularity which is not shared with the other two box 
topologies, and which only appears when the transverse momenta of the final states in 
the centre of momentum frame, Q\ and Q2, are zero. This singularity is known as the 
double parton scattering (DPS) singularity [2], and it is associated with all of the loop 
particles becoming on-shell and collinear with the initial state gluons. The reason why the 
singularity is known as such is that it corresponds to the physical process in which two 
gluons each split to produce an on-shell, collinear quark-antiquark pair, and then the four 
resultant partons interact to produce two Z bosons. The four partons interact in pairs from 
different gluons in two separate annihilation interactions, which is essentially the definition 
of a double parton scattering. 

None of these singularities are restricted to the box diagrams. The conditions to 
have a collinear or soft singularity in a one-loop diagram are well-documented [3, 4]. The 
double parton scattering singularity will occur for any one-loop diagram which satisfies 
the following criteria. First, the two initial state particles must be massless, and each 
of these initial state particles should be connected to two loop particles which are also 
massless. Then, the four massless loop particles should interact in two separate pairs, with 
particles from different initial state particles interacting. There is no restriction on the final 
state from each interaction, only that it should have total invariant mass squared which is 
timelike. Such diagrams will also generically contain collinear singularities. 

The collinear and soft singularities in any loop diagram of the Standard Model are 
completely cancelled by the numerator structure. This has to be the case - otherwise the 
loop diagrams would be infinite for arbitrary values of the external invariants, and one 
would get nonsensical infinite values for the differential cross section of various processes. 
To what extent is the DPS singularity cancelled in Standard Model loops? 

This answer to this question is of interest to two groups of people. The first of these is 
the NLO multileg community, who need to know where the singularities are in a loop inte- 
gral, and how bad they are, to ensure (for example) accurate numerical evaluation of the 
loop integral [2, 5, 6]. It is also of interest to the multiple parton interaction (MPI) com- 
munity, since the nature of the DPS singularity in one loop diagrams determines whether 
part of these diagrams should be regarded as a leading order (LO) double parton scattering 
process or not. If the infra-red DPS singularity in one-loop diagrams of the appropriate 
structure is not integrable at the cross section level, then the singularity should be absorbed 
into multiparton distributions, and a part of each loop diagram can be associated with LO 
DPS. On the other hand, if the DPS singularity is integrable at the cross section level, then 
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the one-loop diagrams can simply be regarded as contributing to single parton scattering 
processes. 

The framework suggested by Snigirev for calculating the cross section for double parton 
scattering processes [7] anticipates that there should be an unintegrable DPS singularity 
in one-loop diagrams of the appropriate structure. We show this explicitly in the following 
few paragraphs. Note that from henceforth, we shall refer to the framework of [7] as the 
'dPDF framework' (for reasons that will become clear shortly). 

Very generally, one may write the cross section for the process pp — > AB + X via 
double parton scattering as follows: 
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The a symbols represent parton-level cross sections. Tij(x\, x 2 , b; Q A , Q 2 B ) is the 
impact-parameter space two-parton GPD (6-space 2pGPD), whilst T ij (xi,x 2 ,r;Q 2 A ,Q 2 B ) 
is the transverse momentum space 2pGPD (fc-space 2pGPD). Tij{x\, x 2 , b; Q A , Q 2 B ) has a 
probability interpretation as the probability to find a pair of quarks in the proton with 
flavours ij, momentum fractions x\x 2 , and separated by impact parameter b, at scales 
Qa and Qb respectively [8, 9]. The fe-space 2pGPD is the Fourier transform of this with 
respect to b, and has no probability interpretation. 

Under the dPDF framework, it is assumed that IYj(xi, x 2 , b; Q A , Q 2 B ) may be approx- 
imately factorised into the product of longitudinal and transverse pieces, with the longitu- 
dinal piece being given by the double PDF (dPDF) object of [10] for the case in which the 
two scales Qa and Qb are equal. The transverse piece is typically taken to be flavour and 
scale independent: 
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Figure 2. A diagram that apparently contributes to the leading order DPS cross section according 
to the dPDF framework (see text). The black circles are hard processes, the grey blobs are protons, 
and the lines are partons. 

The first two terms on the right hand side of (1.3) are associated with changes in the 
dPDF due to independent branching processes - processes in which there are a pair of 
partons, one of which has the appropriate x and flavour, and the other of which splits, 
either giving rise to the other parton of the appropriate x and flavour, or removing it. 
The final inhomogeneous term represents the increase in the dPDF due to a single parton 
with momentum fraction x\ + X2 splitting into a pair with the appropriate x values and 
flavours. We call this the 'sPDF feed term' for obvious reasons. The functions Pj^ f j 1 j 2 (x) 
that appear in this term are known as the 1 — > 2 splitting functions, and may be obtained 
trivially at LO from the real splitting parts of the usual splitting functions. 

Say we wish to calculate the cross section for a DPS process for which Qa = Qb 
(A = W ± , B = W ± , for example). It is clear that if we use (1.2) in (1.1) with the dPDFs 
at arbitrary scale being obtained from some inputs at a low scale according to (1.3), then 
the result for the cross section will contain a term which contains the accumulated sPDF 
feed parts of two dPDFs being multiplied together. Pictorially, this term corresponds to a 
sum of terms with the structure of figure 2. 

Since the dPDF framework includes figure 2 in the LO DPS cross section, this frame- 
work predicts that the loop process of figure 2 should contain a piece which is proportional 
to [as log(Q 2 )] N at the cross section level (where iV is the sum total of branchings that 
occur on either side of the two hard processes A and B). For such a piece to exist, each 
branching in the diagram has to be associated with a transverse momentum integration 
J d\k\/\k\ at the cross section level - even the vertices at which the branches on either side 
of the diagram split into two. The leading log part of the cross section is then associated 
with the region in which the transverse momenta are strictly ordered along the branchings 
on either side of the diagram. 

The (/ d\k\/\k\) 2 factor that is to be assigned to the two 1 — > 2 branching vertices 
in the diagram must originate in the DPS singularity of the loop integral, which must 
therefore be infra-red divergent like a logarithm squared at the cross section level prior to 
subtractions. At the amplitude level, the DPS singularity in the loop integral would have 
to be (log(py)/p^) a5 to give such a divergence at the cross section level (where px is equal 
to the absolute value of the transverse momentum of all of the particles on the upper (or 
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the lower) leg of the loop, in the centre of momentum frame of the two particles initiating 
the loop). 

On the other hand, explicit calculations of certain four- and six-point amplitudes within 
the Standard Model, which contain diagrams of the appropriate character, show that these 
diagrams at least do not contain unintegrable DPS singularities [11-15]. The general 
impression one gets from the literature is that the DPS singularity in any Standard Model 
one-loop diagram is completely cancelled, such that the loop integral is finite in the limit 
Pt -» 0. 

In this paper, we present a detailed and general study of the DPS singularity in one- 
loop integrals. To begin with, we only focus on the four-point diagram that can contain 
a DPS singularity - i.e. the crossed box. In section 2, we present results for the DPS 
singular parts of certain crossed box diagrams, which have been obtained using traditional 
loop integral techniques. Some of these have been extracted from the available literature, 
whilst others are derived by us. We show that in certain Standard Model crossed box 
diagrams, the DPS singularity is not completely cancelled, but is instead relegated to an 
integrable logarithm. 

It is not efficient to evaluate the DPS singular part of a given four-point (or re-point) 
diagram using traditional loop integral techniques, since such techniques involve the full 
calculation of the rest of the integral (which we are not interested in here). Further, we 
can gain little insight from these techniques as to why a particular box integral has a 
DPS divergence of a given nature (in particular, why the Standard Model boxes have a 
DPS divergence which is at most a logarithm of Qz)- In section 3 we derive a framework 
for evaluating the DPS singular part of a crossed box diagram which only requires the 
calculation of simple leading order light-cone wavefunctions and tree-level matrix elements. 
Using this framework, we reproduce and provide physical explanations for all of the crossed 
box results found in section 2. 

The framework that we derive for calculating the DPS divergence of a crossed box 
has the advantage that it is very easily generalised to loops with arbitrary numbers of 
external particles. In section 4, we use the generalised framework to check and generalise 
the results of [14, 15] for the DPS divergence in six-photon helicity amplitudes. We also 
use it to determine the nature of the DPS singularity in figure 2, and compare the result 
with the predictions of the dPDF framework. Based on the outcome of this comparison, we 
comment on the validity of the dPDF framework versus the 'two-parton GPD' framework 
for calculating proton-proton DPS recently introduced by Diehl and Schafer [9]. 

2 Singularities in the Crossed Box 

We consider a generic crossed box diagram with the momenta and helicities labelled as in 
figure 3 (note that any or all of the helicities could be zero in general). For the moment, 
we do not specify the nature of the external and loop particles. We do however impose the 
conditions that are necessary for the crossed box to contain a DPS singularity - namely, 
that the incoming particles (with momenta p\ and P2) should be on-shell and massless, 
whilst the outgoing particles (with momenta Q2 and Q2) should either be on-shell and 
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Figure 3. The crossed box topology, with annotations that demonstrate our labelling conventions 
for the particle names, helicity and momenta. The particle names are written in bold in this figure, 
whilst the helicity labels are accompanied by grey arrows. The arrows on the lines merely indicate 
the direction of momentum flow, and do not necessarily signify a fermion line. The thin lines 
represent masslcss particles, whereas the thick lines represent particles with invariant mass squared 
equal to M 2 . 



massive, or off-shell such that Q\ and Q\ > 0. For the purposes of calculational simplicity, 
the squared four momenta of Q\ and Q2 shall be taken to be equal in all of the boxes 
studied. The common four momentum squared Q\ = Q 2 will be denoted by M 2 . Further, 
we work at all times in the centre of momentum frame, and choose the z axis to be aligned 
with the spatial part of p\. We define: 



s = {pi+ p 2 ) 2 t=(pi- Qi 
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The d dimensional loop integral associated with the crossed box has the following generic 
form: 
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The nature of the external and loop particles determines the numerator factor TV, but 
not the denominator. L is defined such that M only includes the trace structure of the 
crossed box amplitude, and does not include overall factors such as coupling constants and 
colour factors. For future reference, we write here the numerator factors for each of the 
specific crossed boxes that we will consider as examples in this paper, and which are drawn 
in figure 4: 
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Figure 4. The crossed boxes that we shall consider as examples in this paper. The helical lines 
represent gluons, the wavy lines Z bosons, the dashed lines scalars, and the lines with arrows 
represent fermions. 



where e\ is the polarisation vector corresponding to helicity A (A = ±1 for gluons, and ±1, 
for Z bosons). The numerator factor for figure 4(a) written above is of course not the initial 
expression you would write down, which would contain contain factors of v q + a g 7 5 before 
^* and (where we use the notation of [12] - v q (a q ) is the vector (axial) coupling of the 
quarks in the loop to Z bosons). However, the terms in this initial expression proportional 
to Vqa q cannot contribute to the loop integral according to charge conjugation invariance 
[12], whilst the terms proportional to v q and a 2 , can both be shown to have the trace 
structure written above (some anticommutation of the 7 5 matrices plus use of (7 s ) 2 = 1 is 
required in the latter case). Thus, the numerator of the gg — > ZZ loop integral is equal to 
the above trace structure, times some overall coupling constant which we drop. 

We recall that a crossed box contains collinear and double parton scattering Landau 
singularities. In the case of </> 3 theory in four dimensions, these nominal singularities both 
correspond to actual infinite values of the integral. We may calculate the most divergent 
part of the crossed box in this theory (drawn in figure 4(d)) using the elegant method 
presented in [4]. Note that we must work in 4 + 2e dimensions to regulate the collinear 
divergence which appears for arbitrary values of the kinematic invariants. The result of 
the calculation is: 



£</>,4D = ~~ 7~~L 77Z\ — 1~ l ess divergent terms 



4vr 3 
(ut - M 4 )e 

47T 3 

= 5 — h less divergent terms (2-4) 

sQ 2 e 

where Q2 2 is the transverse momentum squared of the second massive particle (= ut—M 4 = 
Qi 2 )- We note that the result obtained here corresponds to minus the expression presented 
in equation (4.101) of [4] - the latter result is out by a factor of (-1) because the two poles in 
equation (4.99) of [4] are on the wrong sides of the real axis. One observes the appearance 
of a factor 1/e in the expression (2.4) which corresponds to the collinear singularity and is 
infinite in the limit e — > 0, and a factor I/Q2 2 which corresponds to the DPS singularity 
and is infinite in the limit Q2 — > 0. A critical point to note is that the DPS singular part 
of the 4D scalar crossed box is not integrable - that is, if ones takes its modulus squared 
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and integrates it over the final state phase space, then one obtains an infinite contribution 
to the cross section (the result is proportional to (IQ2 2 IQi^\ 

In more complex four dimensional theories, there exists the possibility that the collinear 
and DPS singularities in crossed box integrals may exhibit less singular behaviour, due to 
the fact that there is now a nontrivial numerator factor M which may vanish at the singular 
points. Indeed, this seems to be the case for every Standard Model crossed box one can 
construct obeying the appropriate conditions (i.e. p\ = p| = 0, massless particles in the 
loop, and Q\ and Q\ > 0). In all of these boxes, the collinear divergence vanishes, and the 
DPS singularity is relegated to a logarithm in Q2 2 at most. 

Let us give some examples of this logarithmic behaviour of the DPS singularity drawing 
from the established literature. The first example we shall consider is gg — > HH via a 
massless quark loop. Glover and van der Bij have calculated the crossed box integral for this 
process [11]. However, they only present results for the helicity matrix elements calculated 
using a general quark mass m q in the loop. The helicity matrix elements are equal to the 
sum over loop integrals for the six different loop topologies contributing to gg — > HH, 
all multiplied by various factors (coupling constants, colour factors). We can nevertheless 
extract the leading low Q2 behaviour of a single gg — > HH crossed box from these results 
as follows. First, we strip the multiplying factors from the helicity matrix elements to 
obtain expressions for the loop integrals summed over topologies. This turns out to require 
extreme care since the authors of [11] have chosen to factor some constants out of their 
matrix elements and into their expression for da/dt. Then, we take the expansions for 
the scalar loop integrals in the low m q limit (found in Appendix B of [12]), insert them 
into these expressions, and take m q — > to obtain the sums over loop topologies for the 
massless quark case 1 . Such a limit is perfectly well defined since the loop integrals do not 
contain collinear singularities. Finally, we isolate any low Q2 divergences in the resulting 
expressions - these can be equated to twice the leading low Q2 behaviour of the relevant 
crossed boxes. The reason for this is that only the crossed box loop integral can contain 
a DPS singularity, and there are two crossed box topologies that contribute equally to 
99 -> HH. 

Performing this procedure, we find that there are two helicity configurations for which 

the crossed box diverges as Q2 — > - these are the ++ and configurations. The 

corresponding leading low Q2 behaviour of the crossed box integral with either of these 
helicity configurations is: 

W ++ ) = W-) = - 8M > 3 '°^> (2-5) 

s 

This result may be directly checked by decomposing the gg — > HH crossed box loop 
integral to scalar integrals using FeynCalc [17], before inserting the low m q expansions of 
the scalar integrals and taking m q — > 0. We obtain the same expression using this method. 

: An alternative approach would be to insert the dimensionally regulated scalar loop integrals with 
massless internal lines found in [3] . One has to exercise some care in analytically continuing the results of 
[3] to the region of present interest however - see the detailed discussion in [16]. Of course, either method 
gives the same result. 
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Another example of a crossed box satisfying the appropriate conditions is gg — > ZZ, 
again via a massless quark loop. Glover and van der Bij have calculated the loop integrals 
for this process as well [12]. As in the gg — > HH case, they only present amplitudes 
summed over all box topologies and for the case of general quark mass - however we can 
extract the leading low Q2 behaviour of the m q = gg — )• ZZ crossed boxes from these 
results using the same technique as was applied in the gg — > HH case. 

We denote the helicity configuration in a crossed box integral by AiA2//i/i2 where Ai 
and A2 correspond to the helicities of gluons 1 and 2, and /ii and fi2 correspond to the 

helicities of Z bosons 1 and 2. Then only the + + ++, , + H and h+ 

integrals are divergent in the limit Q2 —> 0: 
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Unfortunately we cannot check this result using FeynCalc as it requires the Passarino- 
Veltman reduction [18] of tensor integrals of index 4, which FeynCalc cannot handle. We 
remark in passing that the same results for Lp,ps are obtained if the final state Z bosons 
are replaced by off-shell photons (with Q 2 = Q 2 = M 2 ) or W bosons, except that Mz in 
(2.6) should be replaced by M or M\y. The coupling constant factor that multiplies Lp,ps 
in the full expression for the amplitude is (v 2 + a 2 )g 2 z g 2 in the gg — > ZZ case. The coupling 
constant factor for gg — > 7*7* may be obtained from this by setting v q = Q q ,a q = and 
replacing gz by e, whilst that for gg — > W + W~ is obtained by setting v q = —a q = 1 and 
replacing g z by g w /(2\/2). 

Despite assertions to the contrary that exist in the literature [13], some of the Standard 
Model crossed box loop integrals are divergent in the limit Q2 - > 0. However, a key point is 
that they are not sufficiently divergent to cause the cross section for gg — > ZZ or gg — > HH 
to diverge (J dQ^ 2 log 2 (Q2 2 ) = finite). In particular, the log squared singularity in the 
cross section that would be anticipated by Snigirev [7] does not exist. 

The phenomenon by which the numerator of Standard Model crossed box integrals 
always vanishes at the DPS singular point such that the DPS singularity is integrable 
demands explanation. In [4, 15], it is suggested that 'gauge dynamics' is the physical 
cause. If this is the case, then we should not see similar behaviour in a 4 dimensional theory 
which does not have a gauge symmetry - for example, scalar gluon theory (also known as 
massless Yukawa theory). Therefore, to test the hypothesis of [4], we examined the crossed 
box loop integral associated with the process g s g s — > g*g* in scalar gluon theory, where 
the final state scalars are off-shell by the same timelike amount. The Feynman diagram 
corresponding to the integral is figure 4(c). 

To calculate the leading low Q2 behaviour of the g s g s — > g*g* crossed box loop integral, 
we use two methods. First, we perform the Passarino-Veltman reduction of the integral 
'by hand' in Maple [19], before inserting the expansions of the scalar integrals for small 
loop particle mass found in [12], and then taking the limit of zero loop particle mass. The 
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other approach involves using FeynCalc to perform the Passarino-Veltman reduction. Both 
approaches return the same result: 



Ldps = 4tt 3 log(Q 2 2 ) (2.8) 

We see that the singular behaviour of this box is exactly the same as the Standard 
Model boxes - i.e. the DPS singularity becomes an integrable logarithm (and the collinear 
singularity disappears). This example indicates that we cannot uniquely associate a loga- 
rithmic DPS singularity with gauge theories, and therefore that gauge symmetries are not 
the direct driving force behind this behaviour. 

Although the scalar crossed box integral in four dimensions has a very different singular 
behaviour to the Standard Model and Yukawa boxes, the same integral in six dimensions 
(corresponding to 6D <ft 3 theory) has exactly the same singular behaviour as the 4D SM and 
Yukawa boxes. We can calculate the most singular part of the 6D scalar box by applying 
the method found in section 4.6.2 of [4] to D = 6. In this case, we do not need to deform 
the number of dimensions to D = 6 + 2e since there are no collinear singularities in the 
integral. 

Repeating the steps (4.95)-(4.97) of [4] with D = 6, we obtain: 

3 • f 1 dad(3 

J [sa + u — M z — ie\p + [{u — M z )a — u — ie\ 

where a and (3 are Feynman parameters. It is simple to perform the integration over /3, 
which gives: 

f 1 da 

Ldps^D = J o sa + u _ a M 2_ ie Hsa -M 2 + (u- M 2 )a - ie) (2.10) 

- ln((n - M 2 )a -u- ie)] 

The real part of this integral is finite as Q2 — > 0, and so for the purposes of extracting 
the leading low Q2 singularity we can ignore it: 

f 1 da 

Ldps&qd-^i I — — tttt [— «vrG(Af 2 — sa — (u — M 2 )a) + inQ(u — (u — M 2 )a)l 

J sa + u — M z 1 J 

Jj^sm (2 . u) 

s 

As asserted, the 6D scalar box has a logarithmic DPS singularity in its crossed box. 

There must exist some characteristic that is common to 6D scalar theory, 4D scalar 
gluon theory, and the Standard Model that ensures that the leading DPS singularity in 
these theories is converted from a single inverse power of Q2 2 to a logarithm (and that 
the collinear singularity vanishes). Using traditional techniques for handling loops, it is 
exceedingly difficult to elucidate the mechanism by which this occurs. The reason for 
this is that we lose contact with the original structure of the loop integral when we start 
introducing Feynman parameters (and, in the Yukawa and SM cases, even before this when 
we perform the Passarino-Veltman reduction). In the next section, we shall introduce a 
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technique for directly calculating the portion of a crossed box loop integral which contains 
the DPS singular point when Q2 = (i.e. the point at which all of the internal lines go on 
shell). As the evaluation of the portion of the integral is direct, neither Passarino-Veltman 
reduction nor introduction of Feynman parameters needs to be performed. By use of this 
method, we will discover the physical origin of the logarithmic DPS singularity in 6D scalar 
theory, 4D scalar gluon theory, and the Standard Model. 



3 Physical Investigation of the Crossed Box 

We would like to investigate the nature and origin of the part of the amplitude L which 
is most singular as the transverse momenta of the produced particles goes to zero. This 
part of the amplitude is associated with the region of the loop integration in which the 
transverse part of the loop variable, k, is small (i.e. much less than y/s and masses of 
produced particles). The reason for this is that, when the transverse momentum of the 
produced particles is zero, the small k region contains the point in which all four of the 
loop particles go on shell simultaneously. 

Therefore we study the contribution to L coming from the small k region in the case in 
which the transverse momenta of the produced particles is also small. The method we use 
is similar to that described in section V of [20], although we fix some errors and address 
some subtleties of which the author of [20] did not seem to be aware. 

To begin, we apply a Sudakov decomposition to all of the vectors in (2.2). We define 
lightlike vectors n and p as follows: 

p= -±=(1,0,0,1) n= 0,0,-1) (3.1) 

An arbitary four vector A may be written in terms of these vectors plus a transverse 
part A (which only has x and y components) as follows: 

A = A + p + A~n + A (3.2) 

such that: 

AB = A + B~+A~B + - A B (3.3) 

Writing out all of the four momenta in (2.2) in terms of n, p, and a transverse part, 
(2.2) becomes: 

I d d ~ 2 kdk + dk- 



(2k+k~ - k 2 + ie)[2(k+ - Qi)(k~ - Q2) - (k - Q 2 ) 2 + ie] 

X [2(k+ + Qf)(k~ - Q2 ) - (k - Q2) 2 + ie][2k+(k~ -Qi- Qf) - k 2 + ie] 

In deriving (3.4), we have used the fact that, in our chosen reference frame for which 
pi oc p, P2 oc n, conservation of four momentum implies: 

pi = (Qf + Qt)p P2 = (Or + Q~2> (3-5) 
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In the following discussion, an important point to bear in mind is that Qf and are 
always positive (provided the masses of the produced particles are not zero). 

The part of L that we are interested in is the low k portion of the integral, which we 
shall denote as Lp,ps- Our strategy to evaluate Lp>ps will be to perform the k~ , k + and 
k integrals in that order, making copious use of the fact that k in the integration is small. 

The k~ integration is straightforward. When < k + < Qf, only the second k~ pole 
in the denominator lies on the upper half complex plane, so we close the contour on the 
upper half plane and pick up the pole at: 

+ (3 ' 6) 

When — Qf < k + < 0, only the third k~ pole in the denominator is located in the 
lower half complex plane, so in this case we close the contour on the lower half plane and 
pick up the pole at: 

^« + |^ (3 ' 7) 

Finally, when k + < — Qf or k + > Qf, all of the poles lie on one side of the real axis, 
so we close the contour on the other side and get zero for the value of the integral. Putting 
it all together, we find that the result of the k~ integration is the following: 

f f° N \b--b- 
L DPS = -2iri d d ~ 2 k dk + -, — - s r (3.8) 



(k++Ql) 

1 



2(-Of " Qf)(k - Q2) 2 (-2fc+Qr + fc ff;gj 2 - k 2 + ie 
[ Q * N \k—k- 

+ 27TZ / (fk / dk + 



\k\<Qt,Q7 J o (2k+Qz + k -^E$f ~k 2 + ie 

1 

X 2(Qf + Qf)(k - Q 2 f {-2k+Ql + fc ff:gj 2 - k 2 + ie 

We note that the terms k + (k - Q 2 ) 2 /(k + + Qf) and k + (k - Q 2 ) 2 /(k + - Qf) appear 
in some of the denominator factors. These terms are negligible except where k + ~ — Qf 
or k + = Qf . However, the region of k + which is relevant to the leading Q 2 singularity in 
L is \k + \ <C Qf iQT ■ This is because, when Q 2 vanishes, the configuration in which all of 
the loop particles are on shell corresponds to k = Q^ 71 (i- e - k + = 0). Therefore, for the 
purposes of finding the leading singularity in L, we can drop the k + (k — Q 2 ) 2 / (k + + Qf) 
and k + (k — Q 2 ) 2 / (k + — Qf) terms in the denominator. For similar reasons, we can replace 
and k% by Q^ and set k + = in the numerator. Then, Lp>ps becomes: 

2m f d d ~ 2 k 
LDPS ^2{Qf + Qf) l\ k \ <<Qt , Q - (k - Q 2 ) 2 (3 - 9) 

dk + M \k-=Q~,k+=0 



Q + {2k+Qz - k 2 + ie) (-2k+Q7 - k 2 + ie) 
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Given that \k\ <C Qf,Q^[, the integrand of the k + integration in (3.9) is strongly 
peaked near the origin, and falls off rapidly before either of the two endpoints of integration 
are reached. We can replace the limits of the integration by ±00 without affecting the 
leading singularity in the integral. This allows us to perform the k + integral using contour 
integration, closing in the lower half plane and picking up the pole at k + = k 2 / '(2Q^): 



Ldps 



(^Y f ^- 2 fcA/y =Q -, fc+=0 

'<Qt + Qt)(Qt + Q2) J\k\«Qj,Ql (k - Q 2 ) 2 k 2 



Noticing that 4(Qj + Q 2 ){Q\ + Q 2 ) is simply equal to 2s, we obtain a compact 
expression for the leading Q 2 singularity in L: 

Ldps -— ~ — / 77 — , 2> 2 (3.11) 

2s J\k\«Q+,Q- (k - Q 2 ) 2 k z 

The same result may be obtained by closing the k + integration in the upper half plane. 

Using (3.11), we can reproduce the leading low Q 2 behaviour of all of the DPS boxes 
described in the previous section. To obtain the 4D scalar box result (2.4) , we set N = 1 
and d = 4 + 2e (note that, just as in section 2, we must perform the calculation here in 
slightly more than 4 dimensions to regulate the collinear divergence in the loop integral): 

(2vri) 2 f d 2+2t k 

Ldps^ad =— — / Ti ^ \9i 2 (3.12) 

2s J\k\<.Qi,Qr (k - Q 2 ) 2 k 2 

(2vri) 2 f d 2+2£ k 



2s J {k-Q 2 ) 2 k 2 

4tt 3 
seQ 2 2 

We can expand the domain of integration to infinity because the integrand is strongly 
peaked at k = when Q 2 is small. The usual method of Feynman parameters has been 
used to arrive at the final result. 

The form of the integrand in (3.12) makes particularly clear the interplay between 
the collinear and DPS divergences in the 4D scalar box integral, and their origins. When 
Q 2 7^ 0, there are effectively two distinct poles in the k integration, producing an overall 
logarithmic divergence in the integral. One of these is associated with the loop particles on 
the right hand side of figure 3 becoming collinear (fc = 0) whilst the other is associated with 
the particles on the left hand side becoming collinear (fc — Q 2 = 0). As Q2 is reduced to 
zero, the two poles merge to form a double pole and the divergence in the integral becomes 
stronger (single inverse power rather than logarithmic). The double pole is now associated 
with all of the particles in the loop becoming collinear, and the stronger divergence in the 
integral is precisely the DPS divergence. 

Let us next consider the 6D scalar box: 

(2vri) 2 f d 4 k 

L DPS,<{>,6D =— ^ / —. — , 9> 2 (3.13) 

2s J\kKQ+,Q~ (k - Q 2 ) 2 k 2 
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From this expression, we can clearly see that the 6D scalar box does not possess any 
collinear divergences. 

In the 6D case, we cannot straightforwardly apply the method of extending the in- 
tegration region to infinity that we used in the 4D case. The reason for this is that the 
integrand no longer falls away sufficiently quickly as \k\ —> oo, and we would get infinity if 
we extended the integration region. 

We could evaluate the integral (3.13) by imposing a sharp cutoff A on the integration 
over |fc|, where |Q 2 | « A C Qf,Q~. However, in practical terms it is simpler to use 
dimensional regularisation to extract the leading low Q2 behaviour in Ldps , ,0,6D- We 
evaluate the integral in 6 — 2e dimensions - this allows us to extend the integration region 
to infinity without getting an infinite result. The previously infinite contribution from the 
high k end of the integral now manifests itself as a term containing a single pole in 1/e. 
This can simply be dropped, since we only want the contribution from the low k end of 
the integral. Indeed, we discard every term except for the most singular term in Q 2 . As is 
typical, the dimensional regularisation approach is conceptually more difficult to handle - 
but it produces the same result as the sharp cut-off for the leading singularity in Q 2 . 

Applying the method, we obtain a result which agrees with (2.11): 

{2mf f d*-*k 27r 4 log(Q 2 2 ) 
Ldpsmd =— J (fc _ Qa)2fc2 " UV pole * (3.14) 

It turns out that we must use the dimensional regularisation method to evaluate the 
integral (3.11) for the Yukawa (fi<f> — > <f>*(fi*, SM gg — > HH and SM gg — > ZZ cases as well. 
We evaluated the numerator factors for these integrals using FORM [21]. The results of the 
calculations are listed below - for the gg — > HH and gg — > ZZ cases, we only list results 
for the helicity amplitudes which give a nonzero result for Lops (i- e - are divergent in the 
limit Q 2 -> 0): 

Yukawa g s g s -> g*g*: 

(2m) 2 f d 2 ^kT,[(Q 2 7+ - #)(Qr7 + + #)(QiV ~t + (J 2 )(Q^7 _ + # " 2 )] 

^DPS ~ 



2s J (k- Q 2 ) 2 k 2 

(3.15) 

UV pole 
3i„„/n.2\ 



~47T d log(Q 2 

SM gg -> HH: 



W ++ ) = W-)^- 8M » j3 '°^ 2 > (3.16) 

s 
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SM gg -»• ZZ: 

4vr 

Ldps(+ + ++) = ^DPS = ( ) = — 



s - 2M| + sJl - 4M|/s log(Q 2 



(3.17) 



4vr 3 

^DPs( + + ) = ^DPS = (--+ + ) = 



s - 2M| - s Jl - AM 2 /s log(Q 2 



2\ 



(3.18) 

The results (3.15)-(3.17) agree with those presented in section 2, both in terms of 
dependence on kinematical variables, and in terms of the numerical pref actors. 

In the numerator factor of each of these box integrals, the terms with the smallest 
number of powers of k and/or k — Q2 are proportional to k ■ (k — Q2) - the coefficients 
of the terms with lower powers of k and/or k — Q2 are all zero. A consequence of the 
numerators having this structure is that the leading Q2 singularity in each amplitude is 
demoted from a single inverse power of Q2 to a logarithm. A further consequence is that 
the amplitudes are free from collinear singularities. 

Let us consider the broad features of the method that we have just introduced for iso- 
lating the low Q2 singularity of a box. It consists of performing two sequential integrations 
over the full real axis, picking up the contribution from exactly one pole each time, and 
then performing the integration over k. Picking up the contribution for a particular pole is 
equivalent to replacing the denominator factor corresponding to the pole by a delta func- 
tion (x2-7ri). Essentially, our method is equivalent to replacing the k 2 and (k — Q2) 2 factors 
in the denominator by 2iti5{k 2 ) and 2iri5[(k — Q2) 2 } respectively, before multiplying by — 1. 
We then neglect all the numerator terms during the k integration other than the ones with 
the lowest powers of k and/or k — Q2, such that we pick up the leading singularity in Q-^. 
It is not hard to show the equivalence explicitly: 

- (2vri) 2 / d d k- N0 —. ir , — 5 -5(k 2 )S((k - Q 2 ) 2 ) (3.19) 

V ' J [{Pi + k - Q 2 ) 2 + ie)[{p 2 - k) 2 + ie] y 1 yy ^ ' V ' 



(2vri) 2 / d d - 2 kdk + dk- — 
J [iPi 



+ k — Q 2 ) 2 + ie] [(pa - k) 2 + ie] 



2k~ V 2k ~J 2[(jfe+ -QJ) - (fc- -Qz)k+/k-} V 2(A;+-Q+) 

-\2ki) 2 [ d d ~ 2 kdk + dk- 5-^ 5 ( k+ = —) 5(k- = Qz+ (fc ~ Q2) l \ 

K ' J 2sk 2 (k-Q 2 ) 2 V 2fc-y V 2 2(k+-Q+)J 

+ higher order in Q2 

(2ni) 2 f , d _ 2 , ^k+=o,k-=Q- , ,. , , „ 
' d a fc^r- ; + higher order in Q 2 



2s J k 2 (k - Q 2 ) 2 

It should not be a surprise that the leading Q2 singularity of a box can be obtained 
by replacing the k 2 and (k — Q2) 2 denominator factors by delta functions. Notice that 
the leading Q2 singularity always appears in the real part of L. This corresponds to the 
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Figure 5. The two cuts of the crossed box that can give rise to on-shell particles. 

imaginary part of a box amplitude M since L is always multiplied by — i (along with vertex 
factors etc.) to make an amplitude. But we can obtain the imaginary part of an amplitude 
by using the Cutkosky rules [22, 23]. Thus, twice the real part of L is given by minus the 
sum over all cuts for which the cut propagators may be put on shell (the minus comes from 
the fact that Ai oc —iL). There are two such cuts for the box diagram, which we have 
drawn in figure 5. They give equivalent contributions in the small k and k — Q2 limit, 
with both contributions being equal to minus (3.19). Putting everything together, we see 
that the Cutkosky rules predict that the leading Q2 singularity in the real part of L (= 
the leading singularity in L) is given by (3.19). 

Inserting the values for d and N for the 6D scalar, 4D Yukawa, and 4D Standard Model 
crossed boxes into (3.19), one might get the impression that the real parts of the loop 
integrals L for these boxes are all ultraviolet divergent. It is well-known that an ultraviolet 
divergence exists in the 4D Yukawa and Standard Model crossed boxes - however, this 
occurs in the imaginary part, as can be verified by examining the loop integral expressions 
in the large k limit, and remembering that that a factor of i appears during Wick rotation. 
What we have not written down explicitly in (3.19), but is easy to show, is that for large k 
both delta functions cannot be satisfied simultaneously. Therefore the integral is effectively 
cut off at large k and the real part of L for any crossed box is UV finite. The appropriate 
integration region in k space for the real part of L is an ellipse with the foci at and Q2 
and semi-minor axis length M/2. This approximates to a circle of radius M/2 centred at 
the origin when Q2 is small. 

In the presence of the two delta functions, the remainder of the integrand in (3.19) 
can be decomposed into two factors, corresponding to the two Feynman diagrams of figure 
6(b). Given that the lines with momentum p\ + k — Q2 and P2 — k are almost on shell when 
k and k — Q2 are small, we can use completeness relations to further decompose the upper 
diagram of figure 6(b) into three smaller diagrams (divided by two propagator factors) - 
see figure 6(c). This procedure is very similar to, say, the textbook decomposition of the 
matrix element for e~X —> jY into e~ — > 7e~, e~X — > Y (divided by a propagator factor) 
in the collinear limit (see Chapter 17 of [24]). 

Applying the decomposition of figure 6, the leading low Q2 divergence of a general 
crossed box may be written as follows (recall that our labelling conventions are given in 
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Figure 6. Decomposition of the box integrand. 

figure 3): 

L dps {\iMm2) = E / d d k5(k 2 )5((k - Q 2 ) 2 )^ L S 2 2 L S ^P2;P2 - k,k) (3.20) 

x^lTM-'P i + k-Q 2 ,Q2- k)M s l^B(k, Q2 ~ k; Q 2 ) 
xM2XSa(pi + k-Q 2 ,P2- k; Qt) 

^a~>bc S2 * s essen tially the light-cone wavefunction to find the pair be with helicities sis 2 
inside the particle a with helicity A [25]. Each of these functions in (3.20) is composed from 
three ingredients - the matrix element from the relevant Feynman diagram in figure 6(c), 
the propagator factor nearest to this diagram, and one further factor R. The last factor 
is equal to the square rooted ratio of the collinear momentum fractions of the upper and 
lower outgoing particles in the relevant Feynman diagram. In the spirit of [24], the matrix 
elements in the <3? factors of (3.20) should be evaluated using the following approximate 
expressions for the loop vectors: 

k = Q^n + k; Q 2 -k = Q+p-(k-Q 2 ) (3.21) 
Pi + k - Q2 = Qfp + (k — Q 2 ); P2 - k = Q±n - k 

■M-^l^b * s ^ ne matrix element for the 'hard process' in which the pair L3L4 with 
helicities S3S4 interact to make particle B with helicity fj,2- Given that we are only interested 
in extracting the leading Q 2 singularity of Ldps, it is actually acceptable to evaluate this 
ingredient of (3.20) with all transverse momenta set to zero. 

It should be pointed out that the formula (3.20) only strictly applies when the masses 
of A and B are equal. The reason for this is that to introduce the R factors which are 
a part of the <!> functions into the box integrand, we have used the fact that R for the 
left hand $ is the reciprocal of the R for the right hand Then we can introduce the R 
functions via 1 = RieftRright- This relation only actually holds when Ma = Mb- In the 
more general case in which Ma is not necessarily equal to Mb-, there will be a prefactor 
equal to M B /M A in front of (3.20). 

From our experience of the QCD light cone wavefunction, we can say that $^^ S2 
will in general factorise into two parts, one of which is only dependent on the transverse 
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momentum of b relative to a, and the other of which is only dependent on the collinear 
fraction of the momentum of a that is carried by b. So, for 3> a -»LiL4 for example: 

^TX 4 (pi;pi + k - q 2 , Q2 - k) = x^ix ( Q +%+ ) K ££u( k ~ Q') ( 3 - 22 ) 

The factor X in (3.22) can be interpreted as the square root of the real splitting part 
of a helicity dependent splitting function. In scalar field theory, the function K^^k) 
is simply the 1/k 2 coming from the propagator factor since the splitting matrix element 
is proportional to 1 in this theory. On the other hand, all of the K functions for QCD, 
QED and scalar gluon theory only diverge like 1/k for small k. The reason for this is 
that all of the 1 — > 2 splittings in these theories are forbidden in the absolute collinear 
limit due to nonconservation of J z . This means that the splitting matrix elements must 
all be proportional to k, which goes together with the 1/k 2 from the propagator factor to 
produce a 1/k dependence for K. In QED/QCD, J z is not conserved for the 5/7 — > qq 
collinear splitting because the initial state must have helicity ±1, and the structure of the 
theory forces the quark and antiquark in the final state to have opposite helicities (i.e. 
total J z = 0). In scalar gluon theory, J z is not conserved for the g s — > qq collinear splitting 
because the initial state has helicity 0, and the structure of the theory in this case forces 
the outgoing fermions to both have the same helicity (i.e. total J z = ±1). 

Recall that we can consider A4 as being independent of k, since we can calculate it in 
the limit in which k is zero. Thus, ignoring M and the X parts of $ in (3.20), which will 
only contribute to the prefactor of the leading Q2 divergence in Ljjps, we can write Lops 
schematically as: 

Ldps^Xwm) ~E / d^kK^uik - Q 2 )/^tX 3 (- fc ) ( 3 - 23 ) 

This integral bears a strong resemblance to the integral that one would use to in- 
vestigate what kind of scaling violations a theory has in its parton distributions, which 
schematically is f d d ~ 2 k\K{k)\ 2 (if this integral diverges logarithmically, then the theory 
has logarithmic scaling violations, etc.). It is reasonably clear from this connection that if a 
theory has logarithmic scaling in its parton distributions, then it will also have logarithmic 
DPS singularities in its crossed boxes. The disappearance of the I/Q2 2 singularity in the 
gg —7- ZZ and gg — > HH boxes may be traced to the fact that the K functions in QCD 
have numerator factors proportional to k. This in turn is caused by the vector nature of 
the QCD theory resulting in a vanishing of the matrix elements for collinear splittings, and 
is not directly related to the gauge nature of the theory as was suggested in [4, 15]. 

For the process g — > qq, we present below explicit expressions for the functions $ 
for all possible helicity configurations. Overall numerical prefactors are emitted in these 
expressions - we give only the dependence on the transverse momentum of the quark k 
and the collinear fraction of the gluon's momentum that goes to the quark x: 

$+4+-(x, k) oc x(e+ ■ k)/k 2 ^^{x, k) oc (1 - x)(e+ • k)/k 2 (3.24) 
<&~Z/~q + {x, k) oc x(e~ ■ k)/k 2 $g^~(x, fc) oc (1 - x)(e~ ■ k)/k 2 



-18- 



e+ (e — ) is the transverse part of the polarisation vector with positive (negative) helicity 
along the gluon direction. 

The unpolarised and polarised g —> q splitting functions divided by \k\ are formed from 
appropriate linear combinations of the mod squares of these functions: 

£ l*£ri a * X2 + %~ X)2 « ^ (3-25) 

A,si,S2 

Let us consider the box integral (3.20) for the process gg — > AB with arbitrary final 
states, and ignore the M and X functions in (3.20) which do not depend on k. If the 
two initial state gluons have the same helicity, then in the limit Q2 = this integral looks 
like / d 2 k(e +z • k){e~ z ■ k)/k A [e +z = (1, i) and e~ z = (1, — i)]. This is logarithmically 
divergent. On the other hand, when the gluon helicities are opposite, we get J d 2 k(e^- z ■ 
k)(e^ zz ■ fc)/fc 4 which evaluates to zero. Thus, the gg — > AB crossed box will not contain 
a logarithmic DPS singularity if the initial state gluons have opposite helicities. It is 
important to emphasise that this statement is totally independent of the final states AB. 
Note this general rule is obeyed for the case of the gg — > ZZ and gg — > HH crossed boxes 
- see (3.16) and (3.17). 

The physical explanation for this phenomenon is as follows. The structure of the QCD 
theory forces the qqqq intermediate state in the crossed box process gg — > qqqq — > AB 
(which is essentially real in the collinear limit) to have total J z = in the collinear limit. 
Then, if the initial state gluons have opposing helicities J z = ±2, there is an issue with 
total J z nonconservation aside from local J z nonconservation at each g — )• qq vertex. This 
manifests itself as a further suppression of the gg — > AB box integral numerator in the 
limit k, Q2 — > 0, which makes the integral convergent. 

If the final state particles AB have spin, then there is one further way in which a 
gg — > AB crossed box can become convergent in the limit Q2 = contrary to naive 
expectations. If the helicities of A and B are such that there is no assignment of helicities 
to the internal lines which simultaneously conserves helicity at the g — >■ qq vertices, and 
conserves J z at the qq —> A and qq —> B vertices in the collinear limit, then the crossed 
box integral will not contain a DPS singularity. The extra numerator suppression in the 
limit k,Q2 —> comes from one or both of the factors M in this case. This rule can be 
seen to hold in the case gg — > ZZ. 

We can make some sense of the prefactors in (3.15) - (3.17) in terms of products 
of square roots of helicity dependent splitting functions using our decomposition (3.20). 
Where the factors M are nonzero, they can only be proportional to M regardless of the 
final state. We also find / d 4 k5(k 2 )5((k - Q 2 ) 2 ) oc / d 2 k/M 2 for small k, Q 2 . Taking d = 4 
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Figure 7. Possible configurations of internal helicity for the gg — > AB crossed box in the collinear 
limit, where we have taken Qf > Qi,Q2 > Qt i both initial state gluons have positive helicity, 
and the final state is (a) a pair of Z bosons with positive helicity (b) a pair of Higgs bosons. Note 
that the permitted internal helicity configurations would be the same for both (a) and (b) if the 
helicities of the gluons were both negative instead. 



in Lops'- 

L D ps(Ai Wa) oc £ X^l S : (*) U " v) (3-27) 

where ^Lil^A ^ s si m ply a function that is equal to 1 if J z is conserved in the 'hard process' 
producing final state A, and zero otherwise. Here x is defined to be equal to Q~± /{Q~t + Q2) 
and y = Q2 /{Q\ + Q%)- Since we have taken the masses of A and B equal and work in 
the centre of mass frame, Qf = Q2 , Qi = Q2 > an d y = x. 

It is clear that the second line of (3.27) provides the factor of log(Q2 2 ), whilst the first 
line provides the prefactor that depends on M and s. Without loss of generality, let us 
take Qf > Q^[, Q2 > Q2 ~ i- e - we take A to be the final state particle that travels along 
the +z axis in the collinear limit. 

Consider the + + ++ helicity configuration for the gg — > ZZ box. In the collinear limit 
there is only one possible assignation of helicities to the internal lines which is permitted 
by the structure of the theory and conserves J z at the hard processes. This is presented 
in figure 7(a). In this diagram, the internal lines with positive helicity both have collinear 
momentum fraction of parent equal to x. Bearing this in mind, and using the formulae 
(3.24), we see that the prefactor of the + + ++ crossed box must be proportional to x 2 , 
which in turn is proportional to (s — 2M| + s</l — 4M|/s)/s. 

This process can be repeated for all other processes and helicity configurations. The 

internal helicity configuration is the same for the gg — > ZZ h+ process - looking 

at (3.24) we can then clearly see that the prefactor must be proportional to (1 — x) 2 oc 
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(s — 2M| — syl — 4Mg/s)/s. For the gg — > HH diagram, the two possible arrangements 
of internal helicities are always the same regardless of the gluon helicities (figure 7(b)). 
For both of these arrangements, the internal lines with the same helicity always have 
complementary momentum fractions (this is different from the gg — > ZZ case, in which 
internal lines with the same helicity always have the same momentum fraction) . As a result 
of this, the prefactor for the gg — > HH process is proportional to x(l — x) oc Mfj/ s. Finally, 
the g s — > qq light cone wavefunctions do not contain any dependence on x (a consequence 
of this being that that P q9s does not depend on x [26]), so the prefactor of the g s g s — > g*g* 
crossed box does not contain any dependence on s or M. 

Actually, we can also justify the prefactor for the 6D scalar box (3.14) using the 
framework of (3.20). For the 6D scalar case A4 is now independent of M, so the equivalent 
expression to (3.27) in this case has a prefactor of 1/M 2 , besides having d 2 k replaced by 
d A k and all helicity labels removed. As there are no complications involving spin for 6D 
scalar theory, we can straightforwardly associate the function X^^^x) with the square 
root of the real splitting part of the (f) — > <f> splitting function in 6D c/> 3 theory, which is 
given by P^(x) oc x(l — x) [27]. Putting everything together, we find that the prefactor 
for the crossed box in 6D 3 theory is proportional to x(l — x)/M 2 = 1/s. 

It is not hard to show that equation (3.27) continues to hold even when the masses of 
A and B are not equal, although one has to bear in mind that x is not necessarily equal to 
y in general. It is easy to use this expression (or the scalar 4D/6D equivalent) to generalise 
the results (3.12) - (3.17) to arbitrary masses for A and B. We only write down one of 
these generalisations here - gg —> AB, where A and B are scalars - and leave the others as 
exercises for the reader. The logQ2 2 prefactors of the DPS divergent graphs in this case 

(++ and ) are identical and proportional to x(l — y) + y(l — x) oc (M\ + M^)/s. The 

two terms in this result are associated with the two diagrams of figure 7(b). 

Let us consider the part of the pp — > AB + X cross section associated with two gluons 
splitting almost collinearly into quark and antiquark pairs, and then these pairs interacting 
to form A and B. Suppressing helicity and colour indices: 

(Tpp^gg->AB+x,DPs(s) = J dXdXf g (X)f g (X)a g g^AB,DPs(s = sXX) (3.28) 
v g ^AB,DPs(s) oc-!- J d A qi d\ 2 5(ql - M 2 )5(q 2 2 - M 2 )5^( qi + q 2 - Pl - p 2 ) (3.29) 
X \LDPS,gg->AB\ 2 

By decomposing Ldps according to (3.20) and (3.22), and then making a few substitu- 
tions for the integration variables in (3. 28), (3. 29), one finds that one can bring (3.28) into 
the form of a double parton scattering cross section expressed in terms of the two-parton 
GPDs T of [9]: 

<?pp^gg^AB+X,DPS oc / dxidxid q q^A{s = xiXis)a qq -^ B {s = x 2 x 2 s) (3.30) 

3 i=l 

f ^ T f I \r \ 
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Figure 8. The six-photon loop diagram. 

The T factors in (3.30) are not the full two-parton GPDs - rather they are the con- 
tributions to these objects coming from a g — )• qq splitting, as defined in section 12 of [9] . 
Equation (3.30) is somewhat schematic, in that our full result (and the full expression for 
the DPS cross section in [9]) is actually a sum over terms containing helicity-dependent 
cross sections and helicity-dependent two-parton GPDs that are either both diagonal or 
both off-diagonal in helicity space. The same formula (3.30) is obtained for the close-to- 
collinear part of the pp — > gg — > AB + X cross section if the masses of A and B are not 
equal. 

Thus, it appears that the close-to-collinear part of the pp — )■ gg — > AB+X cross section 
can be described using double parton scattering formulae involving two-parton GPDs. This 
is an interesting result that is not obvious from looking at the basic expression for the box 
integral (2.2). 

The transverse momentum integrals in (3.30) are not divergent at the infra-red end 
(note that there is actually one transverse momentum integration 'hidden' in each 2pGPD 
factor of (3.30)). This is because the DPS singularity in SM crossed boxes is integrable. If 
one extends all of the transverse momentum integrations to infinity, the result is ultraviolet 
divergent - but one might argue that it is not reasonable to do such a thing. We have seen 
that (3.30) corresponds to the contribution to the cross section from the imaginary part 
of the gg — > AB crossed box matrix element, which is UV finite. The r integral and two 
transverse momentum integrals inside the T factors of (3.30) can be mapped via a linear 
change of variables to the Q% integral and two transverse momentum integrals inside the 
Ldps factors of (3.29). All of the transverse momentum integrations (3.29) are naturally 
cut off at finite values, implying that the ones in (3.30) should be also. 

We would suggest that the key property of SM crossed boxes, which is correctly antic- 
ipated by the two-parton GPD framework, is that they give rise to a cross section which 
is infrared finite. Thus there are no infra-red singularities in these objects which need 
to be absorbed into double parton distributions - they should be regarded entirely as a 
contribution to the single parton scattering process. 
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4 DPS singularity in Loops with More Than Four Legs 

An example of a loop diagram with more than four legs which contains a DPS singularity 
is the six-photon amplitude displayed in figure 8. For the diagram to contain a DPS 
singularity, we must take the initial state particles to be the photons with momenta p\ and 
P4, whilst the remaining particles are in the final state. All of the external particles are 
taken to be on-shell (i.e. pf = 0). The DPS singularity occurs when the total transverse 
momentum of photons 3 and 5 (or equivalently 2 and 6) becomes zero. It is associated 
with the point in the loop integration at which k\ and k% become collinear with whilst 
&3 and Jt4 become collinear with p\. 

The first result for a six-photon helicity amplitude, summed over loop topologies, 
was obtained by Mahlon [28] for the MHV helicity configuration. Since then, numerical 
techniques have been developed for performing the loop integration for arbitrary values of 
the external momentum and helicity [2, 5, 6] and analytical expressions for all of the six 
photon helicity amplitudes have been obtained in [29, 30]. In [14, 15], the behaviour of an 
MHV and an NMHV helicity amplitude (in particular the latter) close to a DPS singular 

point is investigated. The helicity configuration in the MHV amplitude is — Y H Y+ , 

whilst that for the NMHV amplitude is Y +- 1- (the ordering of the helicities here 

corresponds to the numbering of the external momenta, and all helicities are defined relative 
to incoming external momenta). Detailed plots are presented in [14, 15] illustrating the 
approach of the NMHV amplitude to the following phase space point satisfying P% = 0: 

p£ = (-33.5,-15.9,-25.0) pt = (-12.5, 15.3, 22.0) (4.1) 

pt = (12.5, -15.3, -0.3) pi = (33.5, 15.9, 3.3) 

The values given above for each photon are the (x, y, z) components of the four momentum 
- the remaining t component is fixed by the on shellness condition. The momenta p\ and 
~p\ are taken to be along the positive and negative z axis respectively. 

The conclusion drawn from the plots is that the NMHV amplitude is finite at the 
DPS singular point, at least for the configuration of external momenta (4.1). It is also 
inferred that the MHV amplitude is finite at the singular point, from the fact that the 
amplitude does not contain any sharp structure when the Nagy-Soper final state momentum 
configuration [2] is rotated around the y axis (with some rotation angles corresponding to 
quite a close approach to the DPS singular point). It is implied in [15] that this behaviour 
is somewhat surprising, given that simple power counting arguments indicate that the 
amplitude should diverge at the DPS singular point as I/-P5: 2 (similar to (2.4)). 

The loop decomposition technique developed in the last section can be very straight- 
forwardly applied to the present situation, to check the results of [14, 15] and investigate 
in more generality the low Ps behaviour of the six-photon loop integral. One thing we 
can say straight away, bearing in mind our experience with the box integrals and recalling 
that QED is a vector theory like QCD, is that the DPS divergence in the loop integral can 
be no worse that the logarithm of Ps- Thus there is certainly no danger of the 27 — > 47 
cross section being infinite. 
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We can actually reproduce the results of [14, 15] without doing any further calculations. 

The DPS singularity in a particular — h H h+ MHV diagram will cancel when we add 

on all other loop topologies which have their DPS singularity in the same place. The 
reason for this is that when we decompose all of these topologies according to (3.20), and 
then add all of the decomposed integrals together, then one finds that one can extract two 
factors from the result which are equal to the full tree-level matrix element for qq — > 77 
(i.e. the sum of both possible Feynman diagrams). This is, of course, not unexpected. For 
the MHV helicity configuration considered the helicity of both final state photons in these 
matrix elements will be the same. But it is well known that the amplitude for a quark and 
an antiquark to produce two photons with the same helicity is zero (see for example [31]). 

So the leading DPS singularity in the — h H h+ MHV diagram goes to zero - i.e. the 

amplitude is convergent. 

The h ++ NMHV amplitude cannot contain a DPS singularity simply because 

the J z of the initial state is not equal to zero. We showed in the last section that a crossed 
box loop integral with two gluons in the initial state does not contain a DPS singularity 
unless the total J z of the gluons is equal to zero. This result obviously generalises to any 
one-loop integral that can potentially contain a DPS singularity, and still applies when 
the initial state gluons are swapped for photons. Thus, the DPS singularity in the NMHV 
amplitude vanishes on a diagram by diagram basis. 

Note that these results serve as a generalisation of the results of [14, 15] to the case 
of arbitrary initial and final state momenta. Aside from using our loop decomposition 
framework to do this, we can also use it to make some interesting statements about the 
singular behaviour of the other NMHV and MHV amplitudes. First, we can say that no 
NMHV six-photon amplitude can ever contain a logarithmic DPS singularity. The reason 
for this is that, however one distributes the helicities, one always ends up either with the 
initial state photons having opposite helicities, or with one of the pairs of the final state 
photons having the same helicity. On the other hand, there are MHV amplitudes that do 
have logarithmic DPS singularities - for example, the H h ++ configuration. 

In section 1, we saw that according to the 'dPDF framework' for calculating DPS cross 
sections, we should be able to associate a (J d\k\/\k\) 2 factor at the cross section level 
with the DPS singularity of figure 2. However, using the loop decomposition technique 
developed in section 3, we can readily see that no such factor may be associated with the 
DPS singularity of the loop. At most, we will be able to associate an integration of the form 
J d 2 klog 2 (k 2 ) with the DPS singularity at the cross section level, which is not infra-red 
divergent. 

We may deduce that diagrams which have the structure of figure 2 should not be 
included in the leading logarithmic DPS cross section, and thus that the dPDF framework 
of section 1 is not sound theoretically. Given our discussion of the crossed box in the 
last section, this result is hardly a surprise. The dPDF framework incorrectly predicts 
the singular behaviour of even the simple gg —> AB crossed box, and we have to use a 
framework in terms of two-parton GPDs (which have an extra transverse momentum or 
impact parameter argument) to get the behaviour right. 

At the heart of the problem with the dPDF framework is the assumption that the 



- 24 - 



two parton GPD may be factorised into a longitudinal and a transverse piece, with the 
longitudinal piece being given by the dPDF object of [10]. In [9], the evolution of two- 
parton GPDs is carefully studied using perturbative QCD, and it is established that the 
'single PDF feed' does not make any contribution to the two-parton GPD for parton pair 
separations b ^ 0. This means that the two-parton GPD evolves very differently for 6 = 
and b ^ 0, which is inconsistent with the factorisation hypothesis. 

Finally, we must point out that even following the breakdown of the dPDF framework 
for describing proton-proton DPS, the dPDF object of [10] (which we also studied and 
provided explicit leading order forms for in [32]) appears to retain validity as the integral 
of the (colour singlet) two-parton GPD over b [9]. This is not probed directly in proton- 
proton collisions, but is probed in proton-heavy nucleus collisions (see [33, 34]). 

5 Conclusions 

In this paper, we have demonstrated that the DPS singular part of any one-loop diagram 
of the appropriate structure may be simply expressed in terms of the transverse momen- 
tum integral of two light cone wavefunctions and two hard matrix elements. An explicit 
derivation of this expression was given for the four-point case, but it is clear that such an 
expression will continue to be applicable for larger numbers of external particles. 

A naive treatment of Standard Model one- loop diagrams initiated by QED/QCD ver- 
tices connected to massless particles indicates that the DPS singularities in these diagrams 
should be of the same strength as those in the corresponding diagrams with scalars - i.e. 
1/Pt 2 , where Pt is the transverse momentum sum of all of the final state particles on one 
of the loop lines extending between the initial state particles. Using our expression for the 
DPS singularity, we have shown that if the theory that determines the vertices at which 
the initial particles attach has logarithmic scalings in its parton distributions (e.g. QED, 
QCD, scalar gluon theory, 6D cj) 3 ), then the DPS singularity in the one-loop diagram can- 
not be stronger than a logarithm of Pt 2 - There is clearly a suppression of the numerator 
in the Standard Model loops at the DPS singular point which causes their DPS singularity 
to go from 1/pt 2 to log(px 2 ). This can be traced to the vector nature of the QED and 
QCD theories forcing J z nonconservation in, and therefore suppression of, the decay of a 
massless particle to a collinear pair of massless particles. 

We exploited our framework to show that an arbitrary one-loop diagram initiated by 
gluons/photons with fermions running around the loop does not contain a DPS singularity 
if the total J z of the initial state is not zero. The physical reason for this is that the total J z 
of the ffff intermediate state in the loop, which becomes real at the DPS singular point, 
is constrained to have J z = at the DPS singular point by the structure of QED/QCD. 
If initial J z ^ there is then an issue of total J z nonconservation (aside from local J z 
nonconservation at each vertex), which suppresses the loop numerator further at the DPS 
singular point and completely removes the DPS singularity. The DPS singularity in a given 
diagram will also disappear if one or both of the hard matrix elements happen to vanish 
in the limit of collinear, on-shell initial state particles. 
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These general principles were applied to explain why the gg — > ZZ and gg — > HH 
box integrals only contain logarithmic DPS divergences for certain configurations of the 
external helicity. In both cases, a necessary condition for the box to have a DPS divergence 
is that the gluons should have the same helicity (ensuring total J z = 0). In the gg — > ZZ 
case, the Z bosons must have the same helicity otherwise there is no configuration of 
internal helicity which ensures J z conservation at both qq — > Z vertices in the collinear 
limit, and the DPS singularity vanishes. It was shown that the prefactors of log(Q2 2 ) in 
the diagrams with DPS divergences could be rationalised as the products of square rooted 
helicity dependent splitting functions (the prefactors of the scalar gluon and 6D cf) 3 boxes 
could also be understood in this way). 

We also applied our general rules to explain why the particular MHV and NMHV six- 
photon amplitudes discussed in [14, 15] contain no DPS divergence. The MHV amplitude 
does not contain a DPS divergence since its hard matrix elements correspond to diagrams 
in which a qq pair go to two photons of the same helicity, which are zero according to 
the MHV rules for QED. There is no DPS divergence in the NMHV amplitude because 
the total J z of the initial state is not zero. We pointed out that no NMHV six-photon 
diagram can ever contain a DPS divergence, whilst there are MHV helicity amplitudes 
that do contain a DPS divergence. 

Finally, we showed that the behaviour of the DPS singular part of a one-loop amplitude 
(i.e. the close-to-collinear part) is inconsistent with the behaviour anticipated by the 'dPDF 
framework' of Snigirev [7] for calculating proton-proton DPS. On the other hand, it is 
entirely consistent with the '2pGPD framework' of Diehl and Schafer [9]. The root of the 
problem in the dPDF framework is the assumption that a 2pGPD may be approximately 
factorised into a longitudinal and a transverse piece, which is not valid. The dPDF objects 
of [10] are the integrals of the b-space 2pGPDs over b. These are not directly accessed in 
proton-proton DPS - however, they seem to be directly accessed in proton-heavy nucleus 
DPS. 
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